Dissipation in Dynamics of a Moving Contact Line 
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The dynamics of the deformations of a moving contact line is studied assuming two different 
dissipation mechanisms. It is shown that the characteristic relaxation time for a deformation of 
wavelength 27r/|fc| of a contact line moving with velocity v is given as r -1 (fc) = c(«)|fc|. The 
velocity dependence of c(v) is shown to drastically depend on the dissipation mechanism: we find 
c(v) = c(v = 0) — 2v for the case when the dynamics is governed by microscopic jumps of single 
molecules at the tip (Blake mechanism), and c(v) — c(v = 0) — 4v when viscous hydrodynamic 
losses inside the moving liquid wedge dominate (de Gennes mechanism). We thus suggest that 
the debated dominant dissipation mechanism can be experimentally determined using relaxation 
measurements similar to the Ondarcuhu-Veyssie experiment [T. Ondarcuhu and M. Veyssie, Nature 
352, 418 (1991)]. 



I. INTRODUCTION 

Spreading of a liquid on a solid surface usually in- 
volves a rather complex dynamical behavior, which is 
determined by a subtle competition between the mutual 
interfacial energetics of the coexisting phases (the solid, 
the liquid, and the corresponding equilibrium vapor), dis- 
sipation processes and geometrical or chemical irregular- 
ities of the solid surface 0. Interestingly, this dynamics 
can be effectively studied in terms of the dynamics of the 
contact line, which is the common borderline between 
the three phases, by correctly taking into account the 
physical processes in the vicinity of it. 

One of the key issues about this dynamics which has 
remained a subject of controversy is dissipation. There 
are two rival theories in the literature, each depicting 
a different physical picture for the dominant dissipation 
mechanism in the dynamics of partial wetting Q. The 
first approach, which is based on the idea of Yarnold and 
Mason || and was later developed into a quantitative 
theory by Blake and coworkers H, emphasizes the role 
of microscopic jumps of single molecules (from the liquid 
into the vapor) in the immediate vicinity of the contact 
line. The other approach, which was developed by de 
Gennes and coworkers |l]||, asserts that for small values 
of contact angle the dissipation is dominated by viscous 
hydrodynamic losses inside the moving liquid wedge. 

For a partially wetting fluid on sufficiently smooth sub- 
strates, a contact line at equilibrium has a well defined 
contact angle 9 e that is determined by the solid-vapor 
75V and the solid- liquid jsl interfacial energies, and 
the liquid surface tension 7 through Young's relation: 
75V — 7s l = 7cos0 e . For a moving contact line, how- 
ever, the value of the so-called dynamic contact angle 9d 
changes as a function of velocity: Bd > 9 e for an advanc- 
ing contact line and 0d < 6 e for a receding one. Since 



the discrepancy between the two dissipation mechanisms 
appears for small contact angles g] , one can expect that 
receding contact lines are in fact very good candidates for 
experimental determination of the dominant mechanism 
in this regime. 

A classic such example corresponds to wetting of a 
plate that is vertically withdrawn from a liquid at a con- 
stant velocity —v, which was first studied by Landau and 
Levich for complete wetting || . In the case of partial wet- 
ting that was studied by de Gennes |||7j , a steady state is 
achieved in which the liquid will partially wet the plate 
with a nonvanishing dynamic contact angle 9d{v) only 
for pull-out velocities less than a certain critical value v c . 
The dynamic contact angle decreases with increasing v, 
until at the critical velocity the system undergoes a dy- 
namical phase transition in which a macroscopic Landau- 
Levich liquid film, formally corresponding to a vanishing 
9d, will remain on the plate. 

Since the onset of leaving a film occurs at small values 
of contact angle, one can imagine that the two different 
dissipation mechanisms would lead to conflicting predic- 
tions about the transition. In particular, in Blake's pic- 
ture the "order parameter" for the transition 9d would 
vanish continuously as v approaches v c , which makes it 
look like a second order phase transition. On the con- 
trary, de Gennes predicts a jump in the order parameter 
from 8 e /y/3 to zero at the transition, which is the signa- 
ture of a first order phase transition 0(5). This drastic 
difference in the predictions of the two theories can pro- 
vide a reliable venue for testing them. However, such 
experiments have so far proven to be inconclusive due to 
the usual difficulties of tuning into a critical point in the 
presence of disorder || . 

Another notable feature of contact lines is their anoma- 
lous elasticity as noticed by Joanny and de Gennes ||. 
For length scales below the capillary length (which is of 
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the order of 3 mm for water at room temperature) , a con- 
tact line deformation of wavevector k, denoted as h(k) in 
Fourier space, will distort the surface of the liquid over 
a distance |fc| _1 . Assuming that the surface deforms in- 
stantaneously in response to the contact line, the elastic 
energy cost for the deformation can be calculated from 
the surface tension energy stored in the distorted area, 
and is thus proportional to \k\, namely M 



E c ] = 



101 



dk 
2^ 



\k\\h{k)f 



(1) 



The anomalous elasticity leads to interesting equilib- 
rium dynamics, corresponding to when the contact line 
is perturbed from its static position, as studied by de 
Gennes |l(J. Balancing and the dissipation, which 
he assumes for small contact angles is dominated by the 
hydrodynamic dissipation in the liquid nearby the con- 
tact line, he finds that each deformation mode relaxes 
to equilibrium with a characteristic inverse decay time 
r~ l {k) = co|fc|, in which cq = 7^/(3^) where rj is the 
viscosity of the liquid and £ is a logarithmic factor of or- 
der unity The relaxation is thus characterized by 
a linear dispersion relation, which implies that a defor- 
mation in the contact line will decay and propagate at a 
constant velocity cq, as opposed to systems with normal 
line tension elasticity, where the decay and the propaga- 
tion are governed by diffusion. This behavior has been 
observed, and the linear dispersion relation has been pre- 
cisely tested, in a very interesting experiment by Ondar- 
cuhu and Veyssie pl| . 

Here we study the dynamics of the deformations of 
a moving contact line for the two different dissipation 
mechanisms. In particular, we focus on the sizeable 
regime where the contact line is moving, i.e. it is away 
from the depinning transition jl2| , but still not too close 
to the onset of leaving a film. We show that in this 
regime the characteristic relaxation time for a fc-mode 
deformation is given as r~ l {k) = c(u)|fc|. The velocity 
dependence of c{v) is shown to drastically depend on the 
dissipation mechanism: we find c(v) — c{v — 0) — 2v in 
Blake's scheme, whereas in de Gennes' picture c(v) might 
be rather well represented by c(v) ~ c(v = 0) — 4w. We 
thus suggest that monitoring the deformation dynam- 
ics in this regime, along the lines of Ondarcuhu- Veyssie 
experiment pl| , can provide a more practical probe for 
experimental determination of the debated dominant dis- 
sipation mechanism. 

The rest of the paper is organized as follows. In Sec. [ij, 
we discuss the two different dissipation mechanisms and 
derive expressions for the corresponding energy dissipa- 
tion rates. These expressions are then used in Sec. Ill to 
derive the force balance, and thus the governing dynam- 
ical equation. The velocity dependence of the dynamic 
contact angle and a characteristic velocity are studied in 
Sees. [V and [v| correspondingly. While Sec. VI dis- 
cusses the effects of sur face disorder, we conclude with 
some discussions in Sec. VII. 



II. DISSIPATION 



Let us assume that the contact line is directed on av- 
erage along the x axis, and is moving in the y direc- 
tion with an average velocity v, which we assume to be 
positive corresponding to receding contact lines, as in 
Fig. § We can describe the position of the contact line 
along the y axis for any given x and t with the function 
y{x,t) = vt + h(x,t). We further assume that the de- 
formation h(x,t) is only a relatively small perturbation. 
We can now try to evaluate the overall dissipation for the 
deforming contact line within the two different scenarios. 




FIG. 1. The schematics of the system. 



A. Blake approach 



The physical process that is involved in causing dissi- 
pation in Blake's picture, i.e. molecular jumps near the 
contact line, is local in nature Therefore, in any small 
neighborhood the amount of dissipation is completely de- 
termined by the local value of the contact line velocity, 
while all the molecular details of the dissipation is en- 
coded in an effective friction coefficient /it -1 . The overall 
dissipation can then be written as 



(2) 



In the limit of relatively small contact angles, which is 
relevant for our receding contact lines, the inverse fric- 
tion coefficient (or mobility) can be calculated as Q 



fcA 3 
k^f 



exp 



W 
kTF 



(3) 



in which W is an activation energy for molecular hopping, 
A is the distance between hopping sites, k is a character- 
istic "attempt" rate, and fceT is the thermal energy. 
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B. de Gennes approach 

We now focus on the contribution of dissipation that 
comes from the viscous losses in the hydrodynamic flows 
inside the liquid wedge . For a slightly deformed con- 
tact line, we assume that the dissipation can be approx- 
imated by the sum of contributions from wedge-shaped 
slices with local contact angles 6(x, i), as shown in Fig. G]. 
This is a reasonable approximation because most of the 
dissipation is taking place in the singular flows near the 
tip of the wedge @J|[l0| . Using the result for the dissipa- 
tion in a perfect wedge, which is based on the lubrication 
approximation |l|,|l3]| , we can calculate the total dissipa- 
tion as fol 



Ph = 4 / dx 



3f 



6(x,t) 



[v + d t h(x,t)Y 



(4) 



in which I = \n{d max / d min ) with d max given by the size 
of the liquid drop and d m i n being a microscopic length 
scale. The inverse dependence on 9 suggests that for suf- 
ficiently small contact angles the hydrodynamic loss is to 
be dominant 101. 



III. FORCE BALANCE AND DYNAMICAL 
EQUATION 

To find the governing dynamical equation in the long 
time limit, we should balance the total friction force ob- 
tained as 6 (Pi + Ph)/5dth(x, t) with the interfacial force 
7 cos 9(x, t) — (jsv — 7sl) = 7[cos 9(x, t) — cos 9 e ] at each 
point along the contact line. Note that in this section we 
are taking both dissipation mechanisms into account. In 
the limit of small contact angles, we find 



ejxj) 



x[v + d t h(x, t)] = 



7 



9(x.tf}. (5) 



To proceed from here, we need to relate the contact an- 
gle 9(x,t) to the contact line profile h(x,t), which can 
be done through solving for the surface profile of the liq- 
uid drop. One can show that the surface profile z(x,y) 
near the contact line can be found as a solution of the 
Laplace equation (&%. + dy)z(x, y) = 0, so as to minimize 
the surface area. The solution that satisfies the boundary 
condition z(x,h(x,t)) = reads || 



z(x,y) 



y 



dk 
2^ 



h(k,t)e ikx - Wy 



(6) 



from which we obtain 

dz(x,y) 



9(x,t) 



dy 

to the leading order 



y=h(x,t) 



dk 

I - / —\k\h(k,t)e lkx 



(7) 



To the zeroth order, Eq.(||) gives the relation between 
the average dynamic contact angle and the velocity as 



( i \ oM-QD 



(8) 



This relation will be used below to study the onset of the 
transition of the moving liquid drop to a Landau-Levich 
film. 

The dynamical equation (Eq.(|J)), which governs the 
dynamics of the deformation field, can now be written in 
the linear approximation as 



d t h(k,t) = -c(v)\k\h(k,t), 



in Fourier space, where 



= ^9 3 d (v) - 3r)(j£v 
° {V> 3 W £ + 9 d (v) ' 



(9) 



(10) 



with 9d(v) to be found by inverting Eq.(g). The corre- 
sponding form of the dynamical equation in real space 
can be found by Fourier transformation as 



dth(x, t) = —c(v) 



dx' h(x',t) 

7T (x — X 1 ) 2 



(11) 



which reflects the non-locality of the dynamics. 

Relaxation of the contact line's shape while it is mov- 
ing thus takes place with the same dispersion relation 
r _1 (fc) ~ |fc| as a contact line at rest, although the corre- 
sponding characteristic velocity c(v) shows a strong de- 
pendence on the contact line velocity v. 



IV. CONTACT ANGLE-VELOCITY RELATION 

There can be two types of experiments on a moving 
contact line depending on how we prepare it. We can 
fix a value for the contact angle that is different from 9 e , 
and let it move with an adjusted velocity when it reaches 
a steady state. This can be achieved, for example, by 
adding or removing some volume of liquid through a sy- 
ringe that is inserted in a liquid drop at equilibrium. On 
the contrary, we can fix the velocity and let the contact 
angle adjust itself in a steady state. This will be the case, 
for example, when a plate is withdrawn vertically from a 
liquid at a constant velocity. 

Depending on which "ensemble" we are using, we will 
have a fixed value for v or 9d, and we should then use 
Eq.(||) (that relates the velocity and the dynamic con- 
tact angle) to determine the conjugate parameter. The 
term "ensemble" is in fact quite appropriate to use here 
because the two (mechanically) conjugate quantities are, 
in fact, velocity and force- which is determined solely by 
the contact angle. What we have is then either a "con- 
stant velocity" or a "constant force" experiment. It is 
interesting to note that in non-equilibrium systems, in 
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general, different ensembles may not necessarily lead to 
the same result @. 

In this work, we are mostly interested in constant ve- 
locity experiments, and thus we will treat v as a fixed and 
given parameter below unless otherwise specified. We 
will examine Eq.(||) in the limiting cases corresponding 
to the two different dissipation mechanisms and compare 
their predictions. 



A. Blake approach 



The behavior in this regime can be extracted from 
Eq.(||) by taking the limit [lt\ <C Od- Inverting the re- 
sulting equation yields 



2v 
cm 



(12) 



in which c;o = HlO 2 . Note that this holds only for 
v < c/o/2, while Od = identically for v > Qo/2. This 
function is plotted in Fig. ||. 

As can be readily seen from Fig. ^ increasing v would 
lead to decreasing values of 6d until at a critical velocity 
vie = Qo/2 it finally vanishes continuously. A vanishing 
contact angle presumably corresponds to formation of a 
liquid film; a so-called Landau-Levich film. The value 
of the dynamic contact angle Od serves as the order pa- 
rameter for this dynamical phase transition, while v is 
the tuning parameter. The continuous vanishing of the 
order parameter makes the phase transition classified as 
second order. As in the general theory of critical phenom- 
ena, a mean-field exponent /3 = 1/2 is characterizing the 
vanishing of the order parameter in terms of the tuning 
parameter. 




FIG. 2. The reduced order parameter T; = (6d/0 e )i as a 
function of the dimensionless velocity Vi = v/cio for Blake 
mechanism (Eq.(|l^)). The dynamical phase transition at 
Vi c = 1/2 is predicted to be of second order in this picture. 



B. de Gennes approach 

In the opposite limit of [ir\ 3> Od, only the hydrody- 
namic contribution survives, and Eq.(@) leads to 



v — i\J 1 — vA 
( i \ x l' c 



1 

V3 



1/3 



(13) 



in which v = 3\/3«/ c M) and Cho — l^e/i^V^)- 1 



in 




FIG. 3. The reduced order parameter Th = (6d/0 e )h as 
a function of the dimensionless velocity Vh = v/cho for de 
Gennes mechanism (Eq.(^)). The dynamical phase transi- 
tion at Vhc = l/(3\/3) ~ 0.192 is predicted to be of first 
order in this picture. 

The above formula, which holds only for v < c^o/(3v / 3) 
has two branches and only the one that recovers 6^(0) = 
e is acceptable as plotted in FigJ|. While at v = 
Cho/(3V3) we find Od = e /v3, we expect to have Od = 
for higher velocities v > Cho/(3v / 3)- Therefore, the order 
parameter Od experiences a finite jump at the transition 
velocity Vh c = c^o/(3\/3), that is the hallmark of a first 
order phase transition. 



V. CHARACTERISTIC VELOCITY 



Using Eq.fllOD, and 0d(v) that we have found in the pre- 
vious section for the two different cases, we can extract 
the u-dependence of the characteristic velocity. 



A. Blake approach 

We can simplify Eq.(pX|) by taking the limit fir) Od, 

as 



x Note that the expression in Eq.(^3|) is real, and the i is 
retained only to keep the appearance of the formula simpler. 
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Cl(9) = C; 



(14) 



Inserting the form of 8d(v) from Eq.(|12|) yields 

ci(v) = c M - 2v. (15) 

Note that in this approach c(v) is strictly linear in v all 
the way, and it vanishes at the transition point, as plotted 
in Fig. | 




FIG. 4. The reduced characteristic velocity C; = c;/c;o as 
a function of the dimensionless velocity Vi = v/cio for Blake 
mechanism (Eq.Jlil)). The slope of the curve is —2 all the 
way to the transition point where the characteristic velocity 
vanishes. 



FIG. 5. The reduced characteristic velocity Ch = Ch/cho 
as a function of the dimensionless velocity Vh = v/cho for dc 
Gennes mechanism (Eq.(|l7])). The slope of the curve is nearly 
—4 until a square-root singularity sets in near the transition 
point where the characteristic velocity vanishes. 

The characteristic velocity can be well approximated 
by the linear expression 



Ch(v) ~ Cho - 4u, 



(18) 



for a wide range of v, except very near c/jq/(3v3) where 
it experiences a square-root singular behavior as 



Ch(v) ~ C h Q 



V2 



V \ / 1 

I +3 ' 



3x73 c h0 J 



c hQ 



-0 



1 



,3\/3 c h c 



3/2 N 



(19) 



It is interesting to note that although both approaches 
predict a sizeable linear regime for c(v), as manifest in 
Eqs.(p^t) and (18), the corresponding slopes are predicted 
differently. 



VI. SURFACE DISORDER 



B. de Gennes approach 

In the opposite limit of /jij Eq.(^) will be sim- 

plified as 



((i s _ c m ( el e d 



(16) 



Putting in Od{v) from Eq.(|T^) leads to 



/ \ ChO 

Ch{v) = 7=3 



i\fl 



1/3 



1/3 



(17) 



One can again check from this equation that c(v) van- 
ishes at the transition. The above equation is plotted in 
Fig. |. 
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In most practical cases, the dynamics of a contact line 
is affected by the defects and heterogeneities in the sub- 
strate, in addition to dissipation and elasticity that we 
have considered so far. If the interfacial energies jsv 
and 73/, are space dependent with the corresponding av- 
erages being jsv and 7sl, a displacement 8y(x,t) of the 
contact line is going to lead to a change in energy as 



SE d 



dxg{x, vt + h(x, t))5y(x, t), 



(20) 



where 

g(x,y) = lsv{x,y) -j S L(x,y) - (^sv - Isl)- (21) 

Incorporating this contribution in the force balance leads 
to an extra force term g(x,vt) on the right hand side of 
Eq.(|), and thus a noise term on the right hand side of 
Eq.(yl) of the form 



r)(x,t) 



3rjfi£ 



g(x,vt), 



(22) 



to the leading order. Note that this is a good approx- 
imation provided we are well away from the depinning 
transition, and the contact line is moving fast enough 

Assuming that the surface disorder has short range 
correlations with a Gaussian distribution described by 
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(g(x,y)) = o, 

(g(x, y)g(x', y')) = g 2 a 2 5(x - x')5{y - y'), 

we can deduce the distribution of the noise as 

{r)(x,t))=0, 
(r)(x, t)r](x', t')) = 2D(v)6(x - x')8(t - t'), 



where 



D(v) = 



1 2 



[J>0d(v) 



2v \9 d (v)+3 W £ 



(23) 



(24) 



(25) 



A,(») = 



2nc 2 



(v/ci )(l - 2v/c w ) ' 



and 



cri(v) 



(v/ciq)' 



(31) 



(32) 



In the presence of the noise, the contact line undergoes 
dynamical fluctuations. These fluctuations can best be 
characterized by the width of the contact line, which is 
defined as 



The above equations are plotted in Figs. ^ and [?]. Fig. ^ 
shows that the width of the contact line is a symmetric 
function of velocity in this picture, while Fig. u\ denotes 
that the order parameter fluctuations decrease monoton- 
ically with velocity. Note also that these fluctuations 
remain finite at the transition point, which is not typical 
of second order phase transitions. 



W 2 (L,t) 



dx(h(x,t) 2 }. 



(26) 



Using Eq.(Q) with the noise term, we can calculate the 
width of the contact line as 



W 2 (L,t) 



D(v) ^ f n/a dk 

TTC(v) 



■n/L 



1 - e 



-2c(v)\k\t 




_ D(v) j \n[c(v)t/a], ^ < t < Jfa, 

nc(v) \ In (L/a), t » ^, [Z '> 

Similarly, we can study the fluctuations in the order 
parameter field 56(x, t) — 6(x, t) — 6 d . Using Eqs.(|^) and 
(§), we find 



o . i 



0.2 



0.3 



0.4 



0.5 



( 2 2 2 \ 



1 _ e -2c(u)|/s|t 



7TC(W) Jn/L 

vD{v)e 2 {v) ( a 2 \ 

2c(v)a 2 \ 2ir 2 c 2 (v)t 2 J' { ' 

for t » a/c(v). 

The magnitude of the fluctuations of the contact line 
width 



FIG. 6. The reduced width A t = A;/ ( M „ g .° I as a func- 
tion of the dimensionless velocity Vi = v/cio for Blake mech- 
anism (Eq.([n])). It is symmetric with respect to Vi = 1/4. 
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and, correspondingly, that of the order parameter 
nD(v)9 2 d (v) 



a(v) = 



2c{v)a 2 



(29) 



(30) 



are thus both velocity dependent. Again, we expect this 
dependence to be different for the two cases. 



FIG. 7. The reduced order parameter fluctuations 

( 2 2g2 \ 
% — - ) as a function of the dimensionless veloc- 
10 J r—j 

ity Vi — v/cio for Blake mechanism (Eq.(B2)). It is a mono- 
tonically decreasing function of velocity. Note the unusual 
feature that the order parameter fluctuations remain finite 
at the transition point, unlike traditional second order phase 
transitions. 



A. Blake approach 

In this case we have fit] <C d , which together with 
Eqs.fl), (H|), (||), (H), and (gj) yield 



B. de Gennes approach 



Taking the opposite limit fir/ ^> 9 d in Eq. (p5|) , together 
with Eqs.tH), @, ©, and ©, we obtain 



G 



and 



g a v e 

6irr/ 2 £ 2 cl 



(-u - iVT^ 2 ) 1/3 + + iVT - ^) 



1/31 



{-v - iVT~ V s ) 1/3 + (-i/ + i VT^2) 



1/3 



(33) 



7T0 



2/)4 



36?7 2 £ 2 c 2 



/)0 



(-1/ - zVT^2) 1/3 + (-1/ + ^^/T — ^) 



1/31 



1/3 



(34) 



The above equations are plotted in Figs. || and |[ Fig. g 
shows that the width of the contact line is not a sym- 
metric function of velocity in this case. Moreover, the 
order parameter fluctuations do not decrease monoton- 
ically with velocity as shown in Fig. ^. Unlike in the 
previous case, these fluctuations diverge at the transi- 
tion point, which is again not typical of first order phase 
transitions. 



VII. DISCUSSION 
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100r 



FIG. 8. The reduced width A h = A h / 



18Tv V 2 l 2 c 2 h0 



a function of the dimensionless velocity Vh = v/cho for dc 
Gennes mechanism (Eq.(|33|)). Note the asymmetry of the 
plot in this case. 
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FIG. 9. The reduced 

„2 fl 4 



order parameter fluctuations 
Bh = o hi ( 36 ^p^ J as a function of the dimensionless ve- 
locity Vh = v/cho for de Gennes mechanism (Eq.(|34[)). In this 
case, it is not a monotonically decreasing function of velocity. 
Note the unusual feature that the order parameter fluctua- 
tions diverge at the transition point, unlike traditional first 
order phase transitions. 



Because of their anomalous elasticity, contact lines re- 
lax to their equilibrium from an initially distorted config- 
uration with a characteristic inverse decay time r~ 1 (k) = 
c(v)\k\ for each fc-mode. The w-dependence of the char- 
acteristic velocity c(v) is shown to crucially depend on 
the dissipation mechanism, and it can thus be used as an 
experimental probe for the dominant dissipation mecha- 
nism. 

A typical experiment for such investigations is direct 
monitoring of the contact line shape during the relax- 
ation process. If the initial distortion of the contact line 
can be made periodic in a controlled way, like in the ex- 
periment of Ondarcuhu and Veyssie JTl| , one can directly 
map out c(v) and hence determine the dissipation mech- 
anism from its ^-dependence. 

Another possibility is to have relaxation from random 
initial distortions, which will be the case when we pull 
out a naturally rough plate from the liquid. Monitoring 
the dynamics of the contact line in this case will pro- 
vide statistical information about the relaxation process, 
from which one can hope to deduce the relevant features 
discussed in Sec. VI. 



We finally note that this linear theory is not sufficient 
for a complete understanding of the Landau-Levich phase 
transition, and it breaks down upon approaching the 
transition point. This breakdown is particularly mani- 
fest in the divergence that we encountered in the width 
of the contact line at the transition point. To be able to 
have a complete description, one should keep the relevant 
nonlinear terms that can be calculated by extending the 
method of this paper, and resort to perturbative renor- 
malization group approaches for the resulting nonlinear 
stochastic equation. We have performed these studies, 
and the corresponding results will appear elsewhere [fL5| . 

In conclusion, we have studied the relaxation dynam- 
ics of the contact lines, and suggested that monitoring 
this dynamics can provide an experimental probe for the 
debated dominant dissipation mechanism. 
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